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Abstract. LOPES is a digitally read out antenna detected in a subsequent analysis of the radio events 

array consisting of 30 calibrated dipole antennas. It is as a change of the phase at these frequencies. Finally, 

' located at the site of the KASCADE-Grande experi- the dispersion of the analog electronics has been 

^ . ment at Forschungszentrum Karlsruhe and measures measured to account for the frequency dependence 

O ' the radio emission of cosmic ray air showers in of the delay. 

G\ ! the frequency band from 40 to 80 MHz. LOPES Keywords: LOPES timing calibration 
is triggered by KASCADE and uses the KASCADE 

reconstruction of the shower axis as an input for the j Introduction 
analysis of the radio pulses. Thereby LOPES works 

as an interferometer when the signal of all antennas The main part of LOPES (LOfar ProtoypE Station) 

is digitally merged to form a beam into the shower consists of 30 digitally read out, absolutely calibrated, 

direction. To be sensitive to the coherence of the radio inverted V-shape dipole antennas Q), El. 0- The an- 

signal, a precise time calibration with an accuracy tennas are co-located with the KASCADE array with a 

in the order of 1 ns is required. baseline of about 200 m and are triggered about twice 

Thus, it is necessary to know the delay of each per minute by the KASCADE-Grande experiment (4), 

antenna which is time and frequency dependent. [5|. Only a few events per day contain a cosmic ray 

Several calibration measurements are performed to air shower radio pulse which is clearly distinguishable 

correct for this delay in the analysis: The group delay from the noise, as the noise floor is quite high inside 

of every antenna is measured regularly (roughly the KASCADE array. Due to a precise time calibration 

once per year) by recording a test pulse which is of LOPES the digitally measured radio data can be 

emitted at a known time. Furthermore, the delay used to form a beam into the shower arrival direction, 

is monitored continuously by the so called phase Furthermore the cross-correlation of the antennas can be 

calibration method: A beacon (a dipole antenna) calculated to be sensitive to the coherence of the radio 

emits continuously two sine waves at 63.5 MHz and signal Q. This way LOPES is a phased array which can 

68.1 MHz. By that a variation of the delay can be be used as a digital interferometer. 
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Fig. 1: Right: LOPES hardware and timing setup. Left: extension used for delay measurements; a pulse generator 
which is connected to a emitting calibration antenna is triggered by the same sync signal as the LOPES DAQ and 
emits a short calibration pulse. Thus this calibration pulse is recorded with a certain delay within a normal event. 
Due to the setup, the delay of the pulse contains a certain offset, which is the same for every measurement and 
therefore can be ignored when determining relative delays between antennas. 



Whenever LOPES is triggered a trace of 2 6 samples 
with a sampling rate of 80 MHz is read out at each 
antenna with the trigger time roughly in the middle of the 
trace. LOPES is operating in the second Nyquist domain. 
Thus the full information of the radio signal between 
40 MHz and 80 MHz is contained in the data and can be 
retrieved by up-sampling (i.e. the correct interpolation 
between the samples) with the zero-padding method 
0. With reasonable computing time, data can be up- 
sampled to a sample spacing below 0.1ns, so that 
the sample spacing does not contribute significantly to 
uncertainties in the timing. 

To obtain a stable timing of the antenna array the ADC 
clock is centrally generated and distributed via cables 
to all DAQ computers (fig. [TJ. With the exception of 
jumps by full clock cycles (see section HVb . the jitter 
of the clock is negligible. Thus the time calibration of 
LOPES is basically reduced to measure the delay of 
each antenna and the subsequent electronics, i.e. the time 
between the arrival of a radio pulse at the antenna and 
its measurement with the DAQ. This delay is different 
for each antenna (mainly due to different cable lengths). 
As for interferometry only the differences of the arrival 
time of the radio signal between the antennas matters, 
the absolute delay is of minor importance. In this paper 
the term delay is therefore meant in a relative sense. 

The following section explains how the delay is 
measured at LOPES. To achieve the necessary precision 
we also take into account second order effects, like 
the frequency dependence of the delay (dispersion) and 



variations of the delay with time. 

II. Delay measurements 

LOPES as a digital interferometer requires the timing 
precision of the radio pulse in each antenna to be much 
smaller than the period of the filter ringing (ss 17 ns), 
when forming a cross-correlation beam into the arrival 
direction of the cosmic ray air shower. Thus the delay 
has to be known and measured with a precision of about 
1 ns, to be sensitive to the coherence of the radio signal. 

At the beginning of LOPES the delay has been 
measured with solar flare events, but now we have 
developed a new method which does not depend on any 
astronomical sources: Simultaneously with the LOPES 
data acquisition we trigger a pulse generator which is 
connected to a calibration antenna at a defined angle 
and distance to each LOPES antenna. The calibration 
antenna emits a pulse with a known delay after the 
trigger. Therefore after repeating the calibration for 
every antenna the pulse should appear in the data of each 
antenna at exactly the same time, if the delay of each 
LOPES antenna (and the connected analog electronics) 
would be the same. 

However, e.g. due to different cable lengths, this is not 
the case, and indeed the calibration pulses are detected 
at different times for each antenna. The relative delay 
between different LOPES antennas can then be obtained 
by measuring the time differences of the detection times 
of the calibration pulses. 

The pulse detection time can be determined by dif- 
ferent methods: On the one hand a Hilbert envelope is 
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Fig. 2: A short test pulse (FWHM » 10 ns) emitted by a function generator and recorded by the LOPES DAQ. 
The crosses indicate the measurement points sampled by the ADC. The lines show the up-sampled signal and the 
Hilbert envelope of the pulse. In the raw data (a) it can be seen that a part of the pulse is delayed by more than 
100 ns due to the dispersion of the filter. After the correction of the dispersion (b), the pulse is more symmetrical, 
the position of the pulse envelope has moved by about 2 ns, its height (field strength) has increased by about 10 % 
and its FWHM (field strength) decreased by about 10%. The effects on height and width are reduced to a few 
percent, when using the sub-band from 43 MHz to 74 MHz. 



calculated and the pulse position is taken as either the 
position of the maximum or the crossing of half height. 
Thereby the position of the maximum shows slightly less 
jitter (RMS about 0.4 ns) when looking at several events 
recorded within a few minutes and the delays obtained 
from the half height crossing and the maximum are 
consistent within errors. On the other hand the position 
can be defined as the minimum or the maximum of 
the up-sampled trace and the jitter for both of them for 
subsequent events is less than the used sampling spacing 
of 0.1ns (after up-sampling). The delays calculated by 
the maximum and the minimum agree well with each 
other. Thus in principle the delay can be determined 
with a precision of better than 0.5 ns. 

Although the delays determined by the two methods 
(pulse position by the trace itself or its envelope) are 
inconsistent by a few nanosecond, this inconsistency can 
be explained at least partially with the dispersion of 
the analog electronics. This difference is reduced to an 
average of about 2 ns when correcting for this dispersion. 
The remaining inconsistency could not be explained so 
far, and is under investigation. Nevertheless, since the 
calculation of the cross-correlation beam is done with 
the (up-sampled) trace and not with the envelope, we 
currently use the delays deferred directly from the trace 
for the analysis of cosmic ray events. 

III. Correction for the dispersion 

The frequency dependence of the group delay of a 
system is called dispersion. For LOPES SIMfl the dis- 
persion has been measured for the antenna, the analog 
electronics (filter) and the connecting cable 0. The 



dispersion of the cable can be considered negligible and 
the largest contribution to the overall dispersion comes 
from the filter. 

Unfortunately, the dispersion of the LOPES V-shape 
antenna is unknown, because it is not easily measurable. 
But the dispersion of the filter has been measured with 
a network analyzer and can be corrected in the analysis. 
This is especially important as the first 10 LOPES 
antennas are connected to a slightly different filter than 
the later 20 antennas. The successfull correction of the 
dispersion in the analysis software has been proven by 
recording test pulses which are emitted by a function 
generator connected to the RML (fig. [TJ. 

As expected, the filter creates a response only to 
the leading and falling edge of pulses. Thus a delta 
pulse with FWHM << filter width _1 (« 25 ns) should 
be seen as a pulse with a width of w 25 ns (FWHM 
of the power). This pulse is linearly distorted by the 
dispersion such that it is partially delayed by more than 
100 ns (fig. 12}. When correcting for the dispersion this 
distortion can be significantly reduced and the FWHM 
of the pulse power is in the order of 30 ns. In addition 
the pulse is shifted by a few ns and the pulse width and 
height (field strength) change by roughly 10%. Thus the 
correction of the dispersion is necessary for precise time 
and amplitude measurements of radio pulses, and future 
radio experiments should aim to correct the dispersion 
of every system in the signal chain. 

IV. Phase calibration 

The delay of each antenna and the read out electronics 
can be assumed to be roughly constant. However, there 



4 



F. SCHRODER FOR LOPES - TIME CALIBRATION OF LOPES 



Reference (Antenna 4) 
Antenna 13, 63.5 MHz 
Antenna 13, 68.1 MHz 



mm 




500 1000 1500 2000 2500 3000 3500 

eventnumber (first 1 events per day, bad days excluded) 

Fig. 3: Example of the phase calibration: The phase dif- 
ferences at the beacon frequencies between one antenna 
and a chosen reference antenna are shown for the first 
10 events of every day for one year of data taking (08 
May 2008 - 07 May 2009, excluding a few days of down 
time and operational problems). 



are variations on the scale of a few ns over time. In 
addition the LOPES clock distribution is not absolutely 
temperature stable, causing the effective delay of single 
antennas to increase or decrease from time to time by 
integer steps of 12.5 ns (one clock cycle). To correct 
for these steps, originally the phase of a TV transmitter 
inside of the LOPES band was monitored. After the shut 
down of this TV transmitter this method has been further 
developed by setting up a dedicated transmitter (beacon) 
inside of the Forschungszentrum Karlsruhe. 

The beacon continously emits two sine waves of 
63.5 MHz and 68.1MHz at -21dBm each, which have 
a very narrow band width (FWHM < 100 Hz). This 
signal is clearly seen above the noise level by all LOPES 
antennas and the phase can be measured at each of the 
two beacon frequencies at each antenna. For each of the 
two frequencies now the following is valid: As a LOPES 
event contains data from every antenna coincidently, a 
stable delay should lead to a constant difference between 
the phases at the same frequency measured by different 
antennas. Thus a variation in the relative delay between 
two antennas can be detected as a variation in the phase 
difference at the same two antennas at one of the beacon 
frequencies. Subsequently this variations can be taken 
into account in the analysis of radio events. 

In a few events the measurement of the phase is 
disturbed by RFI noise. To avoid corrections for 'un- 
real' variations of the delay and to take into account 
ambiguities of the phase if the difference is larger than 
180° ,a consistency check is performed between the 
phase differences at both frequencies. For most of the 
events the results agree well, and this way the timing 
of the antennas can be monitored and significantly 
improved on a event-by-event basis. 

For the given example (fig. [3]), the correlated drift 
which can be seen in both frequencies corresponds to a 



variation in the effective delay of about 1.5 ns over the 
year and can be clearly distinguished from the jitter of 
the phase differences. The jump in the middle is caused 
by a change of the effective delay by 25 ns during one 
day. The outlier at the beginning is one of the few noisy 
events, for which the phase calibration method would 
fail. 

The accuracy of the phase calibration on short time 
scales is thereby determined by the short-term jitter 
(noise) of the phase measurement. This jitter is a func- 
tion of the trace length and the amplitude of the beacon 
signal. Therefore, for a fixed trace length (2 16 samples at 
LOPES) an emission power of the beacon can be chosen 
such that the phase calibration allows a sufficiently 
accurate correction of the timing. In the case of LOPES 
the short term jitter of the phase differences is in the 
order of 0.3 ns. Furthermore, in figure [3] can be seen 
that there is some additional error on longer time scales, 
as changes in the phase differences are not completely 
correlated between both frequencies. But still, the overall 
accuracy is better than the required timing precision of 
about Ins. 

V. Conclusions 

It could be shown that the delay of each LOPES antenna 
and the corresponding electronic chain can be calibrated 
with a precision of better than 0.5 ns. Nevertheless, 
there is a not fully understood difference between the 
measurement of the time of a radio pulse by either 
looking to the trace or its envelope. In addition, it is 
possible to continously monitor the relative delay of 
each antenna with a precision of below 1 ns. Thus long 
and short term variations of the delay in the order of 
a few nanoseconds can be corrected in the subsequent 
data analysis on an event-by-event basis. 

Furthermore, the dispersion of the analog electronics 
has been shown to slightly affect the time, the width 
and the height of measured radio pulses. Therefore this 
is taken into account in the LOPES standard analysis 
pipeline to produce more accurate results. Also future 
radio experiments should consider the dispersion of the 
whole detection system. 

Finally, the different timing calibration methods of 
LOPES show that radio air shower experiments can 
achieve a relative timing accuracy of better than one 
nanosecond and that this accuracy can be monitored and 
maintained continuously over long periods. 
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